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ABSTRACT 

The shape of the hghtcurve during the rising phase of Type I X-ray bursts is deter- 
mined by many factors including the ignition latitude, the accretion rate, and the 
rotation rate of the star. We develop a phenomenological model of the burst rise 
process and show that simple measures of the burst morphology can be robust di- 
agnostics of ignition latitude and burning regime. We apply our results to the large 
sample of bursts from the Low Mass X-ray Binary 4U 1636-536, and find evidence for 
off-equatorial ignition for many of the bursts. We argue that such behaviour may be 
associated with the transition from hydrogen to helium ignition at accretion rates a 
few percent of Eddington. We show that this model can also explain variations in the 
detectability of burst oscillations, and discuss the implications for other burst sources. 

Key words: binaries: general - stars: individual: 4U 1636-536 - stars: neutron- stars: 
rotation - X-rays: bursts - X-rays: stars 



1 INTRODUCTION 



Neutron stars in Low Mass X-ray Binaries (LMXBs) accrete 
matter from their low mass companions via Roche lobe over- 
flow. Nearly half of these systems show Type I X-ray bursts, 
thermonuclear flashes caused by rapid unstable burning of 
. the accumulating hydrogen or helium after it settles and 
' is compressed on the neutron star surface. A typical X-ray 
burst light curve has a rapid rise (less than 10 s), followed 
by a longer decaying tail persisting for seconds to minutes 
as the star's surface cools. 

The basic properties of Type I X-ray bursts can 
be understood in terms of the stability of hydro- 
s gen and helium burning at different accretion rates 
' (iFuiimoto, H anawa Mivai T1 ll98ll : iFushiki Lamb 1 ll987l : 
' Bildsten 111998. 2QQQ:.Narayan k Hevl II2003I: IWooslev et al 



20041 : ICooper Naravan I 120061 : IPeng. Brown Truran 
20071 ). The boundaries between the different burning regimes 



depend on the local accretion rate m (accretion rate per 
unit area). Below a critical accretion rate rhci ~ 1% of the 
Eddington rate mEdd, bursts are triggered by unstable hy- 
drogen burning. Above this level, hydrogen burns stably via 
the hot CNO cycle, and bursts are triggered by unstable he- 
lium burning once a critical column depth is reached. For 
m > rhci but below a second critical rate mc2, expected 
to be a few percent of rriEdd, the hydrogen at ignition depth 
should burn before helium ignition, leading to a pure helium 
flash. At higher accretion rates, m > mc2, hydrogen does 
not burn completely before the burst is triggered, leading to 
various classes of mixed H/He bursts. Bursting activity will 
cease when the accretion rate is sufficiently high, m > rhcs, 



that helium burning is also stable. The precise levels of the 
critical accretion rates mci,mc2, rrics, are not known pre- 
cisely and will vary depending on factors such as the metal- 
licity of the accrete d material and the heat flux from the 
deep crust and core ([Avasli Joss "1982*: 'Fu shiki &: Lamb ' 
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19871 : iBrown. Bildsten k Ru tledge 1998; Nar avan k Hevl 



There is now observational evidence for bursts in all 
three classes. The global accretion rate M can be inferred 
from the X-ray luminosity, and this should (under the as- 
sumption that accretion is spherically symmetric) predict 
rh and hence the burning regime. Furthermore, within a 
given burning regime, burst rate should rise with accretion 
rate. For a large number of sources, however, these expec- 
tations are not borne out. Transition accretion rates can 
differ from the predicted values by up to an order of mag- 
nitude, and for many sources burst rate actually falls a s 
M increases (|Cornelisse et al. II2003I : iGallowav et"ani2007l l . 
Why this should be the case is not yet clear but i t may in- 
volve either additional stable burning processes (Bilds tenl 
11995 ) or non-spherical accretion and slow fuel spread (so 
that local accretion rate rfi differs from that inferred from 
M) ([Bildsten I [2QQ0I ). Alternatively it may wefl be that X- 
ray luminosity is not a good predictor of accretion rate 
(Hasi nger k van der Klis I 19891 ). For a more in-depth dis- 
cussion of these issues and other burst properties we refer 
the re ader to the recent articles by Istrohmaver k Bildsten I 
(|2006l ) and lGallowav et aL IjaooT"). 

In this paper we concentrate on one particular aspect 
of burst physics - the properties of the burst rise. The shape 
and time scale of the rising portion of the X-ray light curve 
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is controlled by several factors: the point at which the burn- 
ing layer ignites, the emission profile of any burning point 
after ignition (set by the nuclear heating and cooling time 
scales); and the speed at which the burning front propa- 
gates across the stellar surface. Ignition is likely to occur at 
a point and spread, rather than occurring across the whole 
star simultaneously, because of the discrepancy between the 
very short b urning!; time s cales and the much longer accretion 
time scale (|Shara lll982l ). Asymmetric initialization is also 
thought necessary to explain the detection during the rise 
of burst oscillations, variations in brightness on the surface 
of th e neutron star that are modul ated by the rapid stellar 
spin ([Strohmaver Bildstenl liooil. 

On a realistic neutron star, ignition location is un- 
likely to be random due to the various factors that break 
spherical symmetry. Accretion flow may not be spherically 
symm etric - it may occur through an equatorial boundary 
layer ('inogamov Sunvaev 1 1 19991 ) or by channeling onto 
the magnetic poles (in which case the star may manifest as 
an X-ray pulsar). Whether either of these issues can affect 
local ignition conditions or lead to preferred ignition loca- 
tions is not clear, since estimates suggest that fuel should 
spread rapidly across the stellar surface between bursts. 
Another factor that will certainly have an impact is stel- 
lar rotation. Most of the neutron stars that show Type I 
X-ray bursts are thought to be rapidly rotating, which re- 
duces the effective gravity at the equator. Centrifugal ef- 
fects, coupled with the deformation of the neutron star 
due to the rapid rotation, combine to reduce the effective 
gravity at the equator as compared to the poles. This re- 
sults in_a^Joca]_a£creti£n_^a^ higher at the equa- 
tor (jSpitkovskv. Levin Ushomirskvl I2OO2I ) . The column 
depth required for ignition is achieved more rapidly, and 
ignition should occur preferentially at this latitude. In fact 
there are accretion rates where off-eq uator ignition is still 
expected. Cooper & Narayan ' ('2007a') considered the situa- 
tion at high accretion rates m ^ rhcs where helium burning 
is on the verge of stability. There will be a range of accre- 
tion rates (the more rapid the rotation, the larger this range) 
where m > rhcs at the equ ator but not at othe r latit udes. 
Although not discussed by 'Coop er Naravan I (|2007al ) , a 
similar region of off-equator ignition is to be expected at 
m ^ rhci (the transition to stable hydrogen burning). 

Once ignition has occurred, the nuclear burning pro- 
cesses determine the emission from a given point. The 
burning layer will expand during a convective phase, and 
there is a delay before radiative processes take over and 
the Ught curve starts to rise. Light curves from a single 
point (ignoring spreading effects) have been generated by 
a number of authors (Taam 1980; Avasli & Joss 198^; 



Fushiki et al. 1119 92: Taam et al. 1993; Woo slev et al. 112004 



Weinberg, Bildsten Schatz 2006). The shape, time scales 



and strength of the single point light curve can vary sub- 
stantially depending on factors such as the burning regime 
and the composition of the accreted material. There are as 
yet no simple analytic models for this process that take into 
account all of the relevant parameters. 

At the accretion rates of relevance to most burst 
sources, th e burning front is expecte d to propa gate by de- 
flagra ti on (iFrvxell Wooslev Ill982bl : [Hanawa Fuiimoto I 
Il984l: iBildstenI Il995ll ra t her than by detonation 
( Frvxefl Wooslev I Il982^ : IZingale et aLl l200lh . In 



the simplest picture, spreading speed is set by the rate 
at which heat is transpo rted across the burning front (by 
convective processes). ISpitkovskv, Levin Ushomirskvl 
(I2OO2I ) have since shown that rapid rotation will also play 
a signiflcant role: interaction between the uplift (vertical 
expansion) of burning material and a strong Coriolis force 
can act to slow spreading. The degree of asymmetry in 
the spread of the burning front is also relevant to the 
detect ability of burst oscillations in the rise, if they are 
caused by a growing hot spot. 

In this paper we attempt the flrst systematic ex- 
amination of how these various factors interact to af- 
fect the shape and time scale of the rising por- 
tion of the light curve. Previous studies in this area 
have focused on small sample s of bursts such as the 
rare multi peaked bursts ( Bhattacharvva Strohmaver I 
2006a fl), or bursts from the accreting millisecond pulsars 
(Bhattacharyya & StrohmayerJ l2006d ). In this study we 
adopt a much broader remit, motivated by the wide variety 
of shapes exhibited by the bursts of the LMXB 4U 1636-536. 
This particular source is an excellent candidate for this type 
of study: there are over 120 bursts from this source in the 
Rossi X-ray Timing Explorer (RXTE) archive, and it lies 
in a binary with relatively well-constrained properties. By 
comparing the burst properties with the results of param- 
eterized simulations, we show that simple measures of the 
burst shape can be profound diagnostics of ignition latitude 
and burning regime. 

The paper is structured as follows. In Section [2] we de- 
flne simple measures of burst morphology and classify the 
bursts of 4U 1636-536 accordingly. Section [3] gives details of 
the parameterized simulations that we carried out to gen- 
erate model light curves, and looks at the effects on shape 
using the same simple measures. In Section 2] we compare 
the results of our simulations with the data, and consider the 
implications in terms of ignition point and burning regime. 
We present a model that can explain our results (subject 
to the assumptions inherent in our simulations) , and extend 
our analysis to other sources to test its feasibility. We con- 
clude in Section [5] 



2 OBSERVATIONAL DATA ANALYSIS 

The bursting LMXB 4U 1636-536 is a persisten t atoll source 
at a distance of ^ 6 kpc (|Gallowav et al. |[2006h in a 3.8 hour 
binary orbit with a low mass blue star, V80 1 Ara. The bi- 
nary inclination i is in the range 36° — 60° (|Casares et al. I 
I2OO6I ). The bursts h ave bee n studied exte nsively with both 
EXOSAT (|Lewin va n Parad iis Taam l ll993) and RXTE 
(|Gallowav et al. 1 120071 ). Burst oscillations at 581 Hz are 
seen in both normal bursts and superbursts, the latter 
being longer bursts triggered by unstable carbon burnin g 
(jStrohmaver et al. \\l998i : IStrohmaver k M arkwardt I I2OO2I I. 
We include in our stud y all of the burs ts covered by the 
RXTE burst catalogue (jCallowav et al. 1 120071) . The cata- 
logue includes over 1000 bursts observed by the RXTE Pro- 
portional Counter Array (PCA) from December 30th 1995 
to October 29th 2005, 123 of which originate from 4U 1636- 
536. 

In order to study the shape and time scales of the bursts' 
rises we started by constructing light curves from the PCA 
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Figure 1. Examples of kinked and multi-peaked structure in the 
burst rise for bursts from 4U 1636-536. The left panel shows 
Burst 18 for this so urce (using the numbering system from 
I Galloway et al. I (|2007l )): this burst shows a kinked rise. The right 
panel shows Burst 40, which has a strongly multi-peaked rise. 

data, using Standard 1 data (0.125 s resolution) and us- 
ing all available energy channels. For most of the bursts 
this results in a smooth and monotonically increasing light 
curve. For a small group of bursts with peak fluxes below 
5 X 10"^ ergs cm-^ s"^ (bursts 11, 67, 71, 76, 78, 87, 88, 
91 , 118 and 122 for this s ource, using the numbering system 
of iGallowav et al. I (|2007f >). the light curve is not smooth, so 
we exclude these bursts from our analysis. Burst 117 is also 
excluded because the rise is truncated. We also exclude for 
the moment bursts that are either kinked or multi-peaked 
when plotted on a coarser (0.5s) time resolution: that is to 
say bursts where the gradient of the lightcurve drops to zero 
or below during the rise. This group includes the strongly 
mu lti-peaked bursts (40, 4 8, 56 and 111) identified as such 
bv IGallowav et al. I tom as well as a group that are ei- 
ther kinked or weakly multi-peaked (bursts 5, 18, 90, 110 
and 113), see Figure [1] We will discuss these bursts in more 
detail in Section 2] 

In analyzing the light curves of the remaining 103 
bursts we considered several different ways of quantifying 
the shape. The measure that we found to be most pow- 
erful, which we call convexity (C), measures the degree to 
which the light curve is concave or convex. The physical sig- 
nificance of this parameter will become clear in Section (3] 
where we find that it can be used to diagnose the latitude at 
which the burst ignites. Convexity C is calculated as follows. 
We define the burst rise to be the interval where the count 
rate rises from 10 to 90 percent of the maximum count rate 
(corrected for the pre-burst persistent emission). We define 
rise time tr as the duration of that interval. In order to 
compare the shape of bursts of different durations and peak 
count rates we then normalize both quantities so that they 
rise from to 10 in dimensionless units. This process is il- 
lustrated in Figures [2] and [3] Taking d as the re- normalized 
count rate in each bin, and Xi as the identity function (shown 
as a diagonal line in Figure [3)) , we define convexity C as 

N 

C = ^(c. -xOAt (1) 

2 = 
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Figure 2. Light curves at 0.125 s resolution for two bursts 
from 4U 1636-536, one convex and one concave. The left panel 
shows Burst 1 and t he right panel Burst 52 (numbers from 
IGallowav et al. I (|2007l )V The dashed lines mark the portion of 
the burst rise where the count rate lies between 10 and 90 per- 
cent of the maximum (corrected for the pre-burst emission). It is 
this portion of the curve that we use to calculate convexity. 




10 10 

Time (normalized) Time (normalized) 



Figure 3. Re- normalized light curves for the two bursts from Fig- 
ure [2] As before, the left panel shows Burst 1, and the right panel 
Burst 52. Convexity C, defined in equation ([T}, is the integrated 
area above or below the diagonal dashed lines, with areas above 
the line (convex bursts) taken to be positive and areas below the 
line (concave bursts) negative. For Burst 1, we find C = 14.8, and 
for Burst 52 C = -20.5. 

where N is the number of re-normalized time bins and At 
is the re-normalized time bin size. Convexity is effectively 
the integrated area of the curve above or below the diagonal 
line - areas above the line being positive and areas below 
negative. For highly convex bursts, where most points lie 
above the diagonal, C is positive, and for highly concave 
bursts it is negative. Our choice of re-normalization means 
that C lies in the range -50 to 50. 

Figure |4] plots the peak flux Fp of the bursts against 
convexity. At high peak fluxes almost all of the bursts have 
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Figure 4. Peak flux for the bursts of 4U 1636-536 (from 
I Galloway et al. I (|2007l )) against convexity. The bursts divide nat- 
urally into two groups based on peak flux (above or below 
50 X 10 ~^ ergs cm"-^ ^~^)- For convenience in later discussions 
we also split the low peak flux bursts into two groups depending 
on whether convexity is positive or negative. Only three bursts 
(those in the upper left quadrant) sit outside this simple group- 
ing. For these we use the presence or absence of radius expansion 
to determine group membership: Group 3 bursts show radius ex- 
pansion while those in Group 1 do not. The lower dotted line 
shows the minimum peak flux for which we calculate convexity 
(5 X 10~^ ergs cm~^ s— 1). 



Figure 5. Normalized persistent flux 7 (as a fraction of the Ed- 
dington rate, taken from Gallowav et al. (2007)) against rise time 
tr for the 4U 1636-536 bursts. Groups 1 and 2 cluster in the range 
7 ^ 0.02 — 0.04. Group 3 bursts are seen at all values of 7 and 
dominate at higher 7. Groups 1 and 2 have longer rise times than 
Group 3 (for which tr ~ 1 s). There are only four Group 1/2 
bursts at 7 > 0.05. Two (partly obscured) have 7 ~ 0.07. The 
Group 2 burst at 7 = 0.10 8 is Burst 16, one of the anomalous PRE 
bursts identifled in Gallo wav et al. ] (|2006l ) (the other. Burst 18, 
is kinked and is not included on this plot). The Group 1 burst at 
the highest 7 is Burst 19: the burst rise shape is unusual, being 
almost kinked. 



C > 0, whereas at low peak fluxes we find equal numbers 
of bursts with both positive and negative convexities. For 
convenience in later discussions, we assign the bursts to three 
broad groups, as indicated in the Figure. Bursts with < 
50 X 10~^ ergs cm~^ s~^ are assigned to Group 1 (C < 0) 
or Group 2 (C > 0), while the bright photospheric radius 
expansion (PRE) bursts with Fp > 50 x 10~^ ergs cm~^ s~^ 
form Group 3. There are no PRE bursts in Group 1, and 
only one in Group 

Figures [5] and [6] show the relationship b etween accretion 
rate a nd burst rise timjfl tr for the bursts. iGallowav et al. I 
give two parameters that can be used to estimate 
the accretion rate at the time of the burst. The first, the 
normalized persistent flux 7, is the ratio of the persistent 
X-ray flux to the Eddington flux. This is thought to be a 
reasonable estimate of the global accretion rate M as a frac- 
tion of the Eddington rate MEdd- The second measure, Sz^ 
measures position in a colour-colour diagram: Sz is thought 
to increase as accretion rate rises. Both measures show the 
same trend. The bursts in Groups 1 and 2 cluster at low ac- 
cretion rates (a few percent of Eddington) and have longer 
rise times than those in Group 3. The Group 3 bursts have 
short rise times and dominate at higher accretion rates. 

Table [1] summarizes the mean properties of the three 

1 Burst 16, identifled bv lGallowav et al. I (|2006l ) as a PRE burst 
with an anomalously low peak flux, due most likely to the pres- 
ence of hydrogen. 

^ Note that we use a different deflnition of rise time to 
IGallowav et al. 1 (120071) . These authors deflne rise time as the time 
taken for the count rate to rise from 25% to 90% of the peak rate. 
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Figu re 6. Colour-colour diagram position Sz (| Gallowav et al. I 
I2OO7I) against rise time. Groups 1 and 2 dominate at lower 5^, 
while Group 3 dominates at higher values. 



groups of bursts. In addition to the parameters previously 
defined, we include burst fluence Eh (the integrated flux dur- 
ing the whole burst) and the burst time scal e Th (the ratio of 
burst fluence to peak flux), using data from IGallowav et al. I 
([2003). Group 3 bursts have higher fluences and shorter time 
scales, whereas properties for Groups 1 and 2 are similar. 
We al so indicate the p ercentage of bursts in each group for 
which IGallowav et al. I (|2007l ) state that oscillations are de- 
tected during the burst rise phase. Oscillations are far more 
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Property Group 1 Group 2 Group 3 

Number of bursts in group 31 30 42 

TR (s) 3.3 ± 1.0 2.8 ± 1.6 1.2 ± 0.3 

Tfo (s) 10.8 ± 4.9 13.7 ± 6.6 7.7 ± 2.2 

(10-6 ergs cm-2) 0.24 ± 0.02 0.33 ± 0.03 0.48 ± 0.02 

7 0.033 ± 0.003 0.033 ± 0.003 0.052 ± 0.003 

S';, 1.91 ± 0.06 1.69 ± 0.08 2.13 ± 0.03 

Fp (10-9 ergs cm-^ s-^) 22.7 ± 1.8 25.5 ± 1.9 62.7 ± 0.8 

% bursts with oscillations in rise 7% 23% 36% 

Table 1. Mean properties for the three groups of bursts. Most parameters are taken from iGallowav et al. 1 (l2007l ) apart from rise time, 
where we use the definition given in this paper. 



prevalent in Group 3 th an in Groups 1 and 2, as previously 
noted for this source bv lMuno et al. 



3 SIMULATIONS 

To understand what might cause the variations in shape and 
time scale seen in the observational data, we developed a 
simple phenomenological model of the burst rise process and 
ran parameterized simulations to generate light curves. Our 
model consists of three main elements: a time-dependent 
temperature profile used to describe the emission from each 
point on the surface after ignition; a velocity model that 
describes the propagation of the burning front across the 
neutron star surface; and a light curve generation routine 
that models the propagation of photons from the neutron 
star surface towards a distant observer. 

To generate light curves we us e the Oblate 
Schwa rzschild (OS) approximation of iMorsink et al. I 
(|2007l ) to model relativistic light-bending, Doppler shifts 
and gravitational redshift. The OS model, which takes 
into account rotation-induced oblateness, is more ap- 
propriate for very rapidly rotating neutron stars than 
the more usual Schwarzschild + Doppler approximation 
([Poutanen Gierlihski I |2003| ). Assuming that burst os- 
cillation frequency is a good measure of stellar spirjfl? 4U 
1636-536 rotates at ^ 580 Hz. The associated rotational 
deformation is a few percent, depending on the assumed 
mass and nuclear equation of state, lead ing to a small 
but n oticeable influence on the light curve (|Morsink et al. I 
We neglect both special relativistic time delay and 
the additional time delay experienced by initially inwardly 
propagating photons, since these delays are m uch smaller 
(- 10"^ s, IPoutanen Beloborodovl (|2006l )) than the 
time bins we consider (^ 0.1 s). We specify stellar mass M 
and equatorial radius Req and then compute the deformed 
spherical surface using the OS model. 

To start the burst we specify an initial small burning 
area and then track the propagation of the burning front 

^ For the two accreting millisecond pulsars that also show 
burst oscillations, the burst oscillatio n frequency is at or very 
close to the known spi n frequency (|Chakrabartv et al. I l2003l : 
IStrohmaver et al~l |2003| ). The situation for the non-pulsing 
LMXBs is less clear (there are some differ ences in burst oscilla- 
tion p roperties compared to the pulsars, see I Watts Strohmaver I 
(|2006l )). However, in all of the suggested models (Section [T} burst 
oscillation frequency lies within a few Hz of the spin frequency. 



across the star. The stellar surface is divided into a grid of 
patches with area ^ 0.1 km^, and we consider a patch to 
be ignited as soon as the burning front reaches the centre 
point of the patch. Once a patch has started burning, we 
need to specify how its emission varies with time. As dis- 
cussed in Section [1] there are various numerically-generated 
single point emission model s, but no simple analytic mod- 
els. In this study we follow iBhattacharvva Strohmaver I 
(|2006al lbl) and assume that the temperature of the burning 
front follows the following profile after ignition: 

T = To + (Ti -To)[l -exp(-t/tir)], t ^ tm 

— Tmexp(-t/tid), t ^ tm (2) 

where tm is the time at which the temperature reaches its 
maximum Tm = To + 0.99(Ti — To). The time scale tir sets 
the time scale at which the temperature increases, while tia 
sets the time scale on which it decays. Unburnt patches are 
assumed to have a temperature To until ignition. The pa- 
rameters in this model depend primarily on the composition 
of the burning material, which will vary with on accretion 
rate (see Section [Tj). Bursts which are helium-rich, for exam- 
ple, would be expected to have shorter timescales than those 
which contain a higher fraction of hydrogen. Our parameter 
space must therefore be wide enough to take into account 
the expected level of variation. One follow-on question is 
whether this exponential temperature model remains valid 
across all burning regimes. It has been used successfully in 
detailed spectral modelling of bursts at differen t accretion 
rates by IBhattacharvva Strohmaver I (|2006al lbh, but ul- 
timately one would like to see this confirmed by detailed 
nuclear physics calculations. 

We then assume black body emission at the specified 
temperature from each patch. In Figure [71 we show the typ- 
ical single patch light curve . The rise portion is similar in 
shape to the bursts shown in lWooslev et al. 1 (|2004h . The de- 
cay portio n does differ from tha t seen in some of the cases 
studied bv lWooslev et al.1 (|2004h , but because we are focus- 
ing on the rise we never reach the points late in the decay 
phase where the difference would become relevant. This rel- 
atively simple analytic formulation is therefore a reasonable 
approximation to a more realistic calculation, for the short 
rising phase that we are studying. To this we can then apply 
various different beaming functions. 

In adopting this emission model we have made two as- 
sumptions. We have assumed that each point on the star 
will have the same emission profile. This may not be the 
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Figure 7. A typical lightcurve for a single patch on the NS surface 
before taking into account propagation effects. Note that the y- 
axis is linear. The rise time scale ti^ = 0.1 s, and the decay time 
tld = 6 s. All single point lightcurves for the assumed temperature 
profile are intrinsically convex: C ~ 10 for ti^ 1 s, falling as ti^ 
gets shorter. 



case, since emission could vary (as a result, say, of inhomo- 
geneities in fuel deposition or composition). We have also 
neglected any delay between ignition and emission. If a con- 
vective zone develops there will be a delay as the zone grows 
before radiative processes become a more efficient form of 
transpo rt, a nd radia tion starts to escape from the photo- 
sphere ([Weinberg, Bildsten Sz Schatz 2006). If this delay is 
position-independent then there will be no effect on the over- 
all shape of the lightcurve; there could however be some 
position-dependence here that would contribute to varia- 
tions in emission profile across the star. 

In modeling the spreading of the burning front, we use 
th e expression for burning front speed, t 'flame, developed 
bvlSpitkovs kv, Levin Ushomirskv I (|2002[ ). This speed de- 
pends on various nuclear burning parameters as well as the 
(poorly known) strength of frictional coupling between the 
top and bottom of the burning ocean. In general the velocity 
is latitude dependent (fas ter at the equator), because of the 
role o f the Coriolis force. 'Spitkovs kv, Levin Ushomirskv I 
(|2002l ) show that the burning front speed varies as 



1/tfr + ry/tn 



1/2 



(3) 



In this expression, g is the acceleration due to gravity at the 
surface, /ihot is the scale height of the hot burnt material, and 
/ = 20. cos is the Coriolis parameter, O being the angular 
frequency the star and the latitude. The parameter tfr is 
the time scale for frictional coupling between the top and 
bottom layers of the burning ocean. The parameter tn is the 
nuclear time scale of the thermonuclear burning (set by the 
composition of the burning material) , and 77 is a constant of 
order unity. In the case of weak frictional coupling (tfr ^ 

in,l//j3 



COS^ 



(4) 



where the velocity at the pole, Vp ^ y/ghh^/2Qtn- If fric- 
tional effects are stronger (tfr ^ tn), behaviour depends on 
how tfr compares to 1//. Maximum fiame speed is reached 
when tfr = 1//. In this case 



(5) 



with Vpm ^ghhot/^Otn. In the high frictional coupling 
hmit (tfr ^1//) the latitude dependence disappears, and 



Vfla 



(6) 



The shape of the light curve depends on a whole host 
of parameters: however some have a larger effect than oth- 
ers. We start by presenting results for a baseline scenario, 
where some of the parameters are fixed. Later in this sec- 
tion we will show that varying these additional parameters 
has only a limited effect on our findings. In the baseline 
scenario we assume a neutron star mass M = IAMq and 
an equatorial radius i?eq = 12 km. We set the spin rate 
of the neutron star to 580 Hz, and assume a binary incli- 
nation i — 5 0° , in the middle of the range inferred for 4U 
1636-536 bv ICasares et al. 1 (|2006l ) . In the spreading speed 
model we assume that frictional coupling is weak (equation 
lU, so that fiame speed is inversely proportional to cos^. 
We start our simulations by setting an initial burning area 
with r adius 1 km, similar to the expected widt h of the fiame 
front (jSpitkovskv, Levin Ushomirskv ll2002l V In the tem- 
perature model (equation [2]) we fix the following parame- 
ter s;_7o_^J_^eV^j_7j^^_2^8_^ey,a^ = 6 s (as used 
bv iBhattacharvva Strohmaver 1 (|2006al ll. We also assume 
beaming oc cos-^, where is the angle from the normal to 
the surface. 

The parameters that we vary are the ignition latitude 
ttign, the polar velocity Vp and the temperature rise time 
scale tir- We vary ignition latitude in 30° steps from 0° 
(north pole) to 180° (south pole). The parameters Vp and 
tir both depend on the composition of the burning material, 
which will vary with accretion rate as the fraction of hydro- 
gen changes. We choose ranges for these parameters wide 
enough to encompass the level of variation that might be ex- 
pected as we move from pure helium bursts to hydrogen-rich 
bursts (|Wooslev et al. I I2OO4I : [Weinberg. Bildsten Schatz I 
2006). We consider Vp = (1,2,4,8,16) x 10^ cm/s (for ve- 
locities outside this range our simulations give rise times are 
either much shorter or much longer than those observed) . In 
the temperature profile we consider tir = 0.001, 0.1, 0.5, 1, 2 
s, a range wide enough to cover the values generated 
in the single point light curve m odels of, for example, 
[Weinberg. Bildsten Schatz I (|2006l ). 

Figure [8] shows the results from our simulations. It is 
clear that ignition latitude has a dramatic effect on both 
convexity and rise times. We find the following key result: 



Note that this approximation breaks down in a very small re- 



gion near the equator, where 1// - 
use the full expression for velocity. 



In this case we have to 
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Figure 8. Convexity C against rise time tr for simulated bursts, showing the effect of changing Vp and tir- The range of Vp and t\j. 
studied covers the variation expected as the composition of the burning material changes from pure helium to a mix containing a 
substantial fraction of hydrogen. Each panel shows a different ignition latitude: ajgn = 0° (north pole), ajgn = 30°, oiign = 90° (equator), 
ck^ign = 150°, and aign = 180° (south pole). The cases aign = 60° and 120° are not shown but are very similar to the equatorial case. For 
ignition at the equator, convexity is always positive. Convexity decreases as ignition moves towards the poles, and can become negative. 
Negative convexity bursts are much more common for ignition in the northern hemisphere than the southern. Rise times are shortest 
near the equator, increasing as ignition point moves towards the poles, with north pole ignition bursts showing longer rise times than 
south pole ignition bursts. The difference in behavior between the poles reflects the fact that we have set i = 50°. 



bursts ignited on the equator always have C > 0, for the full 
range of Vp and tir considered. Convexity decreases as igni- 
tion latitude moves towards the poles, and can eventually 
become negative at high enough latitudes; negative convex- 
ity being much more easily achieved when ignition starts in 
the northern hemisphere than in the southern. Rise times do 
vary slightly with ignition latitude but depend more strongly 
on spreading speed and the temperature rise timescale. The 
differences between northern and southern hemisphere igni- 
tion are due to the fact that we are observing the star from 
a northern hemisphere vantage point {i = 50°). For bursts 
ignited at the south pole we do not observe the initial slow, 
concave, part of the rise that is visible to us if the burst ig- 
nites at the north pole. A south pole ignited burst therefore 
appears to have a higher convexity, purely by virtue of our 
viewing angle. 

Let us now consider whether our results are robust when 
we vary the parameters that were held constant in our pre- 
vious simulations: neutron star mass and radius, inclination 



angle, and the size of the initial burning area. We also vary 
the remaining temperature parameters, which are functions 
of composition and hence accretion rate. To do this we took 
a representative sample of bursts (Table [2]), varied the other 
parameters, and compared the results to the baseline sce- 
nario with the same ignition latitude, Vp and tir- 

We considered neutron star masses in the range M = 1 
to I.8M0 and equatorial radii i?eq 10 to 14 km, in line 
with reasonable estimates f or a range of equations of state 
(jLattimer Prakash Il2007h . We varied inclination from i = 
0° to 90° (a range much larger than that inferred for 4U 
1636-536, but which will allow us to extend our results to 
other sources in Section |4]) . The radius of the initial burn- 
ing region was varied from 200 m to 1 km. We took t\d in 
the range 4-15 s, and tested the effect of increasing both To 
and Ti by a factor of 5. We also included the effect of the 
RXTE spacecraft response on our simulated bolometric light 
curves. To do this we generated typical PCA response func- 
tions using some of the observed bursts analyzed in Section 
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Burst ID 


«ign n 


Vp (xlO^ cm/s) 


tlr (s) 


1 





2 


0.1 


2 





2 


0.5 


3 





4 


0.1 


4 





4 


0.5 


5 


90 


2 


0.1 


6 


90 


2 


0.5 


7 


90 


8 


0.1 


8 


90 


8 


0.5 


9 


180 


2 


0.1 


10 


180 


2 


0.5 


11 


180 


4 


0.1 


12 


180 


4 


0.5 



Table 2. Fixed parameters for the twelve bursts that we used in 
testing the effects of the parameters that were held constant in the 
baseline scenario. For the equatorial ignition bursts we considered 
a wider range of Vp because these values were a better fit for the 




observed rise times. 



Figure 10. As for Figure [9] but showing the effect on rise time. 
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Figure 9. Convexity for the 12 test bursts (Table [2]). The black 
line shows the baseline scenario. The dashed lines show the max- 
imum deviation from this scenario when we vary ti^? ^5 ig~ 
nition patch size and include the RXTE response function, for 
inclinations in the range 40° — 60° (the range appropriate for 4U 
1636-536). 



[21 We then used XSPEC to fold our simulated light curves 
through the response functions, and re-computed convexity 
and rise times for the folded light curves. 

Figures [9] and [10] show the deviations in (respectively) 
convexity, and burst rise time, as compared to the standard 
scenario, for inclinations in the range 40° — 60° (the range 
appropriate for 4U 1636-536). The results of our previous 
analysis are clearly robust. In Figure [TT] we show results for 
a much wider range of inclinations. The range of variation is 
much wider, but again our main conclusions stand: negative 
convexity requires off-equatorial ignition. 

We also checked the effects of several other parame- 
ters that are not shown in the Figures. One thing that we 
tested was the effect of changing the beaming factor. We 
first tried removing the beaming factor: this had very little 
effect on bursts ignited in the northern hemisphere or at the 
equator. Burst ignited in the southern hemisphere, however. 



Figure 11. As for Figure [9] for inclinations in the range 0° — 
90°. At i = 90°, as expected, north and south polar ignition 
bursts become identical. As inclination falls, convexity increases: 
at i = 0°, for example, the values are too high to explain the 
observations from 4U 1636-536. 



were more prone to having negative convexities. We also 
tried the more physica lly motivated grey atmosphere model 
(|Miller Lamb lllQQSh . In this model, the temperature of 
each patch 



3^4 



■ coa^ + 0.7 



(7) 



where Tgff is the effective temperature of the emitting layer 
and ijj is the angle from the normal. The effect on convexity 
was extremely small compared to the baseline scenario (as 
might be expected, since both have the same -0 dependence). 

We also varied the latitude-dependence of spread- 
ing speed. In the baseline scenario we assumed that we 
were in the low frictional coupling regime, where latitude- 
dependence is strongest (vfiame oc 1/cos^). By the inter- 
mediate friction regime, inflame oc l/cos^^^^, while in the 
highest friction regime the dependence on latitude vanishes. 
We therefore considered two additional models, i;flame = 
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Vp/ cos"*^/^ 6 and Vflame = Vp- Interestingly both of these mod- 
els failed to reproduce the observed convexities. The con- 
stant velocity case produced no negative convexity bursts at 
all. The intermediate friction case did generate some bursts 
with C < but the values were not sufficiently negative to 
match the observations. 



4 DISCUSSION 

4.1 Ignition latitude and burning regime 

We start by summarizing the most important results from 
the preceding sections. In our analysis of the bursts of 4U 
1636-536, we identified three populations: one group with 
higher peak fluxes and two groups with lower peak fluxes. 
The high peak flux group dominate at higher accretion rates 
and have short rise times. All but two of the members of this 
group have positive convexity. The lower peak flux bursts 
divide into two groups of similar size: one with positive con- 
vexity and one with negative convexity. Rise times in this 
group are longer than for the high peak flux group, suggest- 
ing a different burning regime. 

In our simulations, both composition of the burning ma- 
terial (hydrogen fraction) and ignition latitude have a strong 
influence on convexity (Figure [8|. However only bursts that 
ignite near the poles have negative convexities: bursts ig- 
nited on or near the equator always have positive convexity. 
In addition, bursts ignited near the north pole always have 
lower convexity than those ignited near the south pole. We 
can therefore draw the following conclusiorjfl- The Group 1 
bursts, which have negative convexity, must be triggered by 
polar ignition. In the rest of the Section we will outline a 
simple model that might explain why this should be the 
case. 

As discussed in Section [T] there may be a region of 
polar ignition at high accretion rates at the point where 
He ign ition is transitioning to stability (jCooper Naravan I 
l2007al ). However, the same arguments should also apply at 
much lower accretion rates (m ~ rhci), where H ignition is 
on the verge of stability. We therefore propose the following 
model. At the lowest accretion rates, the bursts of 4U 1636- 
536 are mixed H/He bursts, triggered by unstable H ignition. 
The star is however close to the transition where H burning 
stabilizes, in the range of accretion rates where off-equatorial 
ignition is preferred. We expect this range to be reasonably 
broad for 4U 1636-536 because of the rapid stellar rotation. 
As accretion rate increases H ignition stabilize^, and there 

^ Subject, of course, to the assumptions in our modeling. Nega- 
tive convexity bursts could, for example, be generated by equa- 
torial ignition if the single point lightcurves were concave rather 
than convex (see Figure 0. Although we are not aware of any 
physical motivation for such an assumption, we ran several sim- 
ulations with concave single point lightcurves to study the pos- 
sibility. In these simulations, however, we could no longer gener- 
ate positive convexity bursts with long rise times. To explain the 
presence of both positive and negative convexity bursts at low 
accretion rates we would therefore require both types of single 
point lightcurve (temperature profile) to operate simultaneously. 
^ The transition is actually rather more complex. Between the 
regime where H ignition triggers mixed H/He bursts, and the 
higher accretion rate regime where H burning is stable, there ex- 



is a transition to He-ignited bursts. Above the transition 
ignition moves back to the equator. 

Within this picture the Group 1 and Group 2 bursts 
are mixed H/He bursts triggered by off-equatorial H ig- 
nition. We would expect similar numbers of bursts to be 
triggered in the northern and southern hemispheres. We 
therefore suggest that most of the Group 1 bursts, which 
have negative convexities, ignite in the northern hemisphere, 
while most of the Group 2 bursts ignite in the southern 
hemisphere (for which positive convexities are more likely 
The two groups should have similar peak fluxes and rise 
times, in accordance with our observations. The properties 
of the bursts in these groups are in accordance with those 
expected for mixed H/He bursts triggered by H ignition: 
low peak fluxes, rise times of a few seconds, and durations 
longer than 10s. The value of 7 at which the transition takes 
place (7 ^ 0.03) is perhaps a li ttle higher than expected (in 
'Peng. Br own Truran"! (|2007l l. for example, H ignition no 
longer triggers mixed H/He bursts above ~ 1% Eddington), 
but the precise values of the transitions will depend on fac- 
tors such as heating from the deep crust, and 7 is not a 
perfect measure of accretion rate. 

The Group 3 bursts would in this picture be triggered 
by He ignition on the equator. At accretion rates immedi- 
ately above the transition one would expect the bursts to 
be nearly pure He, with the amount of H involved in the 
bursts increasing as accretion rate rises. The properties of 
the Group 3 bursts are in line with those expected for He- 
dominated bursts: high peak fluxes (with radius expansion), 
rise times ^ 1 s, and durations < 10 s. We note that a tran- 
sition to short rise time radius expansion bursts at a few 
perce nt of the Eddington rate was also seen in the EXOSAT 
data fiLewin et al. Ill987h . 

A transition from mixed H/He burning triggered by H 
ignition, to He burning, should result in a drop in burst rate 
and an increase in alpha (the ratio of the energy released 
by stable burning between the bursts to the energy released 
in bursts). If our interpretation is correct this should occur 
when the accretion rate is ^ 3 — 5 % MEdd, and this is 
indeed what is observed (see Figure 16 of iGallowav et alTI 
This picture might also explain some unu sual fea- 
tures of the PRE bursts for this source. Ga llowav et al. I 
(2006) found that while most of the PRE bursts reached 
the Eddington limit for pure He, two had lower peak fluxes, 
requiring some H in the mix. These two exceptional PRE 
bursts occur, as discussed in Section [2l at the highest ac- 
cretion rates. At these rates pure He bursts are no longer 
likely and bursts should have some mixed H/He character 
again, reduci ng the peak flux reached by t he PRE bursts. 
As shown by iGallowav Gumming I (|2006l l , there can be a 



ists a range of accretion rates where hydrogen burns unstably via 
weak f lashes (|Peng, Brown Truran I l2007l : ICooper Naravan I 
"2007b"). The flashes do not trigger associated He burning, so are 
faint and most likely undetectable above the persistent accretion 
luminosity. The weak hydrogen flashes may also show a move 
from equatorial to polar ignition as accretion rate rises. 
^ Some of the Group 2 bursts could ignite on the equator, as this 
would also give rise to positive convexity. This might be expected 
at the lowest accretion rates, and could explain the apparent ex- 
cess of Group 2 bursts at the lowest values of Sz, see Figure [6] 
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substantial percentage of H in the mix before the Eddington 
hmit starts to fall below that expected for pure He. 

The values of convexity that we measured for the bursts 
of 4U 1636-536 ranged from -20 for the Group 1 bursts up to 
+20 for the Group 3 bursts. Although our simulations gen- 
erated bursts with negative convexities in the right range, 
convexities > 10 were harder to generate (Figure [8|). One 
reason for this is that our simulations did not take into ac- 
count PRE (which ah of the Group 3 bursts show). PRE 
tends to flatten the top of the light curve: when we include 
this in our simulations it results in an increase in convexity 
of sufficient magnitude to explain the discrepancy. Changing 
the parameters from the baseline scenario can also increase 
convexity (compare Figures [8] and (9]). A rise in peak temper- 
ature, for example (as might be expected for He-rich bursts) 
increases convexity. 

By comparing the observed values of convexity and rise 
time with those generated by the simulations, we can in- 
fer the range of tir and Vp required to explain the observa- 
tions (comparing Figures |4]|5] and Figure [8|) . For the Group 
3 bursts, for example, we need t\r < 0.5 s independent of Vp, 
while for the Group 2 (and hence also the Group 1) bursts 
we require tir > 0.1 s. We can therefore ask whether the val- 
ues that we infer are in line with the values predicted for the 
suggested burning re gimes. The singl e point lightcurve mod- 
els of .Weinberg. Bildsten &: Schatz I (|2006h predict a tem- 
perature rise timescale of ~ 0.01 s for He bursts, rising to 
^ 0.1 s as H fr action increases. The timescales derived by 
IWooslev et al. ] (2004) using multizone models are slightly 
longer. The limits that we derive are therefore broadly com- 
patible with these models. What is harder to check is the 
validity of the inferred values of Vp. Spreading speed (equa- 
tion [3|) depends on the nuclear timescale, the scale height of 
the burnt atmosphere and the strength of frictional coupUng 
in the burning layers. Naively one would expect the nuclear 
timescale to be related to the temperature rise timescale (al- 
though the degree to which convection develops could skew 
this relationship). The dependence of the frictional coupling 
on the burning regime, however, is very poorly understood. 
Without a better understanding of the burning and spread- 
ing process, it is difficult to say whether the values of Vp sug- 
gested by our simulations accord with the values expected 
for the different burning regimes. 

4.2 Accretion disk effects 

One issue that we have not considered in our modelling 
is the role of the accretion disk. An optically thick accre- 
tion disk extending down to the stellar surface could ob- 
scure the southern hemisphere of the star. At the lowest 
accretion rates, the disk is likely to be truncated sufficiently 
far from the star for ou r unobscured models to be valid 
(jPone Gierlihski lliooi ). As accretion rate increases, how- 
ever, the inner disk is expected to move in towards the stellar 
surface, obscuring the southern hemisphere. 

Radiation pressure from a bright burst may be able 
to push th e disk back, revealing the southern hemisphere 
once more (jShaposhnikov. Titarchuk Habern[2 003) . Even 
if the southern hemisphere were obscured, however, there 
would be little impact on our results. Bursts ignited in the 
northern hemisphere would still have negative convexities, 
and bursts ignited at the equator would still be positive, al- 



though peak countrates would be lower and apparent rise 
times shorter than in our simulations. The largest impact 
would be on bursts ignited in the southern hemisphere: these 
would only be visible to us once the burning front traversed 
the equator, so would all have positive convexities. 

4.3 Burst oscillations 

In Table [1] we summa rized the detectability of burst oscilla- 
tions, as reported bv lGallowav et al~l (|2007l l , for the rising 
phase of bursts from 4U 1636-536. The detectability criteria 
used were based on power exceeding a certain threshold in 
short time bins (0.25s). Oscillations were most likely to be 
detected in Group 3 and then in Group 2. Group 1 bursts 
were far less likely to show oscillations. Can this be explained 
by the pattern of flame spreading if the bursts in Groups 1 
and 2 result from polar ignition, while those in Group 3 re- 
sult from equatorial ignition? Detectability will depend on 
relative amplitude of the asymmetry, the length of time for 
which it persists, and the overall countrate. A long-lived 
lower amplitude asymmetry may be more detectable than 
a short-lived high amplitude asymmetry if the time window 
used for the power spectral analysis is long. 

We carried out a number of simulations to test whether 
our model is compatible with the observations. In Figure [12] 
we show dynamical power spectra from simulated bursts ig- 
nited at different latitudes. To mimic the differences between 
the different groups we assumed a higher peak countrate for 
the equatorial ignition burst. We find that oscillations are 
most detectable in the bright equatorial ignition bursts (our 
model for the Group 3 bursts). For the bursts with lower 
peak countrates, oscillations are more detectable for bursts 
that ignite in the southern hemisphere. This fits with our 
interpretation in which Group 2 bursts are ignited near the 
south pole and Group 1 bursts near the north pole. 

4.4 Multi-peaked and kinked bursts 

4U 1636-536, like many other neutron stars, shows both 
multi-peaked bursts and bursts with kinks (at least two 
points of inflection) in the rise. This group includes those 
that we excluded from our analysis in Section [2] as well as 
bursts with weaker kinks such as Burst 19 (see Figure [5]) and 
the two bursts with negative convexities in Group 3 (Figure 
Sj) . The simulations described in Section O which involved 
single point ignition and subsequent smooth spread, pro- 
duced no simulated light curves with more than one point 
of inflection. 

We therefore extended our simulations to investigate 
the effect of near simultaneous ignition at equivalent lati- 
tudes in the northern and southern hemispheres. We wanted 
to see whether a second burning front igniting before it could 
be engulfed by a burning front spreading from the other 
hemisphere could give a multi-peaked or kinked light curve. 
We ran a series of simulations where north and south pole 
were ignited either simultaneously or at intervals from 0.5 to 
4s. In no case were we able to replicate either multi-peaked 
or kinked shapes. 

We conclude that some additional physics is still re- 
quired to explain the multi- peaked and kinked bursts. 
iBhattacharvva Strohmaver I (|2006al E) put forward a sce- 
nario in which multi-peaked bursts were caused by polar 
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Figure 12. Dynamical power spectra from simulated bursts with 
ignition starting just off the north pole (aign = 10°), at the equa- 
tor (ojign = 90°) and just off the south pole (aign = 170°). Note 
that we start just off the poles because bursts ignited exactly 
at the pole spread in a totally axisymmetric fashion and there 
would be no burst oscillations. For these simulations we used the 
baseline scenario with Vp = 4000 m/s and tij. = 0.1 s. The peak 
countrate for the polar bursts is 2000 cts/s while the peak coun- 
trate for the equatorial burst is 4000 cts/s (see the explanation 
in the text). We use 0.2 5s time windows ( to ma tch those used for 
oscillation detection bv iGallowav et al~l (|2007h ). overlapping by 
0.05s. The contours start at a power of 5 and rise in increments 
of 5. Note the variation in strength and duration. 

ignition in one hemisphere, some kind of stalling on the 
equator, followed by the burning front re- activati ng and 
pro2agating on towards the other pole. In Watts & Maurer I 
(|2007l ) we argued against this model on the grounds that po- 
lar ignition was only expect ed at the highest accretion rates 
(ICooper Naravan |[2007al ). whereas the multi-peak bursts 
are actually seen at low and intermediate rates. However, 
as we have pointed out in this paper, polar ignition is also 
likely to play a role at the transition to stable H burning. 
A mechanism involving polar ignition may therefore still be 
responsible, provided that a viable stalling mechanism can 
be identified. Interestingly three of the multi-peaked bursts 
from 4U 1636-536 sit at the accretion rates where we are pos- 
tulating polar ignition. The fourth, which looks rather differ- 
ent, h owever, sits at slightly higher rates (|Watts &: Maurer"! 
l2007l ). 

4.5 Application to other sources 

If ignition does indeed move towards the pole as we reach 
the upper limits of the accretion rate where H ignition is 
still feasible - and then back to the equator when He igni- 
tion begins - then we would expect to see similar trends in 
other sources. Assuming a favorable inclination we would 
also expect to see positive convexity bursts at the very low- 
est accretion rates (when H ignition should occur at the 
equator), and negative convexity bursts at the very highest 
accretion rates (when He ignition moves off-equator towards 
the pole) . Figure [13] illustrates the type of behavior that we 
expect, and the regions where polar ignition should occur. 
In what follows we will refer to this plot as the F diagram. 
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Figure 13. A schematic of the type of behaviour we expect in a 
plot of 7 (a proxy for accretion rate) against burst rise time for a 
rapidly rotating neutron star: in the text we refer to this as the F 
diagram. At the lowest accretion rates, bursts are triggered by H 
ignition and ignite at the equator. As accretion rate rises ignition 
moves off-equator and rise times increase, forming the lower hor- 
izontal stroke of the F. When H burning stabilises, ignition (now 
by He) moves back to the equator and we move up the vertical 
stroke of the F. Rise times should increase as the proportion of 
H in the bursts increases, and we may move onto the upper hori- 
zontal stroke of the F before He ignition moves off-equator at the 
highest accretion rates. 

We selected a number of additional sources that are 
thought to accrete a mix of H /He from the RXTE burst 
catalogue fCallowav et al. iboOTl ) . We chose the sources that 
span the widest range of accretion rates (those for which 
IGallowav et al. I (|2007l l compute a colour-colour diagram po- 
sition SzjB- We also included two additional sources: the 
eclipsing system EXO 0748-676, which has a large sample 
of bursts at low accretion rates that are thought to be trig- 
gered by H ignition; and the accreting millisecond pulsar 
XTE J1814-338, which has a large burst sample but where 
magnetic confinement may also play a role in setting igni- 
tion latitude. Our simulations can be very easily extended to 
neutron stars with different rotation rates: spreading speed 
is inversely proportional to rotation rate, so a more slowly 
rotating star will have a higher Vp than a star with more 
rapid rotation, all other factors being identical. The change 
in oblateness has a very small effect on convexity and rise 
time. 

We need to be somewhat cautious in comparing differ- 
ent sources, as burst properties will differ. These differences 
will reflect variations in the composition of accreted mate- 
rial, deep crustal heating, and accretion history. The F di- 
agram for one source may be offset from the F diagram for 
another, or different regions of the diagram may be popu- 

^ We excluded 4U 1728-34 because this source is thought to ac- 
crete from an H-poor donor, and XTE 2123-058 because there are 
only six bursts in the catalogue. 



© 0000 RAS, MNRAS 000, 000-000 



12 Maurer & Watts 



lated, depending on the source. A star with an H-poor donor, 
for example, would not trace out the lower portions of the F 
diagram. However, if our model is correct, then there should 
be two transitional accretion rates for each source, for which 
ignition occurs at the pole rather than the equator. As our 
simulations show, not all polar ignition will result in nega- 
tive convexity bursts. However, all negative convexity bursts 
should indicate polar ignition. We therefore expect to find 
negative convexity bursts on the horizontal strokes of the F 
diagram. 

In Figure [14] we plot F diagrams for the other sources 
examined, marking those bursts with negative convexity. As 
for 4U 1636-536 we excluded both the faintest bursts and 
those with kinked or multi-peaked rises. For the accreting 
millisecond pulsar XTE J1814-338, the Eddington luminos- 
ity has not been measured so we show unsealed persistent 
flux rather than 7. The properties of the sources examined 
are summarized in Table O Inclination is for most sources 
unconstrained (although it is expected to be < 70° due to 
the absence of dips or eclipses) . The exception is EXO 0748- 
676, which is a high inclination eclipsing system. We expect 
the range of accretion rates for which polar ignition is im- 
portant to be larger for the more rapidly rotating sources 
(although whether or not we see negative convexity bursts 
would of course depend on the inclination). 

The picture that emerges is certainly not clear-cut. All 
but one of the sources (4U 1702-429) show some negative 
convexity bursts, and many of these sit on what we might 
interpret as the lower horizontal stroke of the F diagram. 
4U 1746-37 is interesting in that it seems to trace out the 
entire upper portion of the F diagram, with two negative 
convexity bursts at high accretion rates that might be polar 
ignition bursts triggered by He burning. There are however 
a few negative convexity bursts at intermediate accretion 
rates (as there were for 4U 1636-536). EXO 0748-676 has 
a large percentage of negative convexity bursts: if these are 
caused by polar ignition then this may seem a little surpris- 
ing, since this source is thought to have a much slower spin 
rate. This source does span a very narrow of accretion rates, 
however, and as a high inclination system we also expect 
the south pole ignition bursts to have negative convexities. 
However, if the accretion rate is not extremely finely- tuned, 
then there may be some other factor (such a strong magnetic 
field) that is controlling fuel deposition and hence ignition 
latitude. The relatively slow spin inferred for this source 
suggests that magnetic field effects might be important. In 
this respect the results from the pulsar (XTE J1814-338) 
are also extremely interesting. This source too shows a large 
number of negative convexity bursts, as one might expect 
if ignition were occurring off-equator at the magnetic pole. 
The phase-locking of persistent pulsations and burst oscilla- 
tions in this source also sug:g:ests that this may be the case 
(jStrohmaver et al. lliooih . 



5 SUMMARY AND CONCLUSIONS 

We have shown, using parameterized simulations, that burst 
rise shape can be a valuable diagnostic of the burning pro- 
cess. Changes in ignition latitude, in particular, can have a 
major impact on burst morphology, and such changes may 
explain the variation in burst rise shape seen in the well- 



studied source 4U 1636-536. We have argued that a change 
from off-equatorial to equatorial ignition might be the hall- 
mark of the transition from H triggered mixed bursts to He 
triggered bursts at low accretion rates. Such a model can 
also plausibly explain variations in the detectability of burst 
oscillations. There are also areas, however, where additional 
physics is clearly required. Our spreading and burning mod- 
els were not, for example, able to generate multi-peaked or 
kinked burst rises. 

Our goal with this work was to develop a simple phe- 
nomenological model to study the interactions between the 
various processes operating during the burst rise, and their 
influence on light curve shape. What we have done is clearly 
simplistic: the various elements of the model will be more 
closely connected than we have considered here, and we have 
made a number of assumptions that may not be valid. The 
parameter space that we have considered, however, is ex- 
tremely wide, giving us confidence in our conclusions. What 
we have done also demonstrates the power of simple mea- 
sures of the burst shape: if our model is valid we have 
been able to identify ignition latitude and rule out latitude- 
independent flame spreading speeds, for example (see also 
"Bhattacharyya & StrohmayerJ (|2007l )l. This type of study 
could and should be repeated as more detailed models of 
the nuclear burning, spreading and emission process become 
available. 
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Source 


Frequency (Hz) 


Bursts analyzed (total in catalogue) 


Number with C < 


4U 1608-52 
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3 
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38 (44) 





4U 1705-44 
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KS 1731-26 
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24 (27) 


1 


4U 1746-37 




20 (30) 


7 


EXO 0748-676 


45 [4] 


72 (83) 


50 


Aql X-1 


549 [5] 


35 (40) 


4 


XTE J1814-338 


314 [6] 


28 (28) 


18 
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[3] ISmith. Morgan Bradtl (Il997h . [4] IVillarreal Strohmaver! (|2004l ). [5] IZhang et al. I (Il998h . [6] IStrohmaver et al. I (|2003l ) 
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